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The pinwheel tiling is the paradigm for a substitution tiling with circular symme- 
try, in the sense that the corresponding autocorrelation is circularly symmetric. 
As a consequence, its diffraction measure is also circularly symmetric, so the pin- 



Q^' wheel diffraction consists of sharp rings and, possibly, a continuous component 



with circular symmetry. We consider some combinatorial properties of the tiles 
and their orientations, and a numerical approach to the diffraction of weighted 
pinwheel point sets. 



■ 1 Introduction 

O ^ 

'. Diffraction is an essential tool for the determination of the atomic structure of solids. The 

discovery of aperiodically ordered materials in the early 1980s flQlllOj added a new dimension 
of complexity to the inverse problem of structure determination [3l H] , and led to a renewed 
, interest in mathematical diffraction theory, which focuses on kinematic diffraction O [6] . It 

also posed the question of what type of ordered structures are possible, and how one would 
be able to detect them. Regular model sets, which are structures derived from a (periodic) 
lattice in a higher-dimensional space, show pure point diffraction; see [llj and references 
^ I therein. However, there are many interesting systems that are not model sets, such as tilings 

with circular symmetry; we refer to [7] for general background on the theory of aperiodic 
order. 

Pinwheel tilings of the plane, first introduced by Conway and Radin [HI [151 [El [13) are 
arguably the simplest examples of substitution tilings with circular symmetry. They are built 
entirely of triangular tiles of edge lengths 1, 2 and \/5- The corresponding right triangles 
occur in both possible orientations, or chiralities, so in fact a pinwheel tiling comprises two 
tile shapes which are mirror images of each other. 

A pinwheel tiling of the plane can be constructed as the fixed point of the inflation rule of 
Figure [H It consists of an inflation with linear scaling factor \/5 and a rotation by an angle 
9? = — arctan(^) ~ —26.565°, followed by a dissection into five triangles of the original size 
and shape, two of the same and three of the opposite chirality. The rotation ensures that 
the original triangle reappears in the inflated patch, and the sequence of iterates therefore 
converges to a space-filling tiling of the plane. At the same time, because (p is not a rational 
multiple of vr, it produces an irrational rotation, which means that in each inflation step one 
obtains triangles in a new direction. In the limit, the set of available directions becomes 
dense on the unit circle, and this produces the circular symmetry of the structure; a proof 
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Figure 1: Inflation rule of the pinwheel tiling. Control points are denoted by dots, with the 
open dot specifying the origin. The inflation rule is compatible with reflection, in the sense 
that the triangle of different chirality is replaced by the mirror image of the patch shown here. 



of the circular symmetry of the autocorrelation can be found in |16l I12j : see also [9] for 
generalisations. 

2 Pinwheel control point sets 

Also shown in Figured] are control points, which form a point set A that is mutually locally 
derivable (MLD) with the pinwheel tiling, in the sense of [8]. The control point of a triangle 
is located at position ^) with respect to the right angle corner of the triangle. In other 
words, the triangle with control point (0, 0) and its short edge along the horizontal axis has 
vertices (— ^), (^,— ^) and (— ^,|), or (— i,— i), (|,— ^) and (^,|), depending on its 
chirality. Given a control point x G ~ C, a triangle can be specified by a triple {x,uj,x), 
where x £ C denotes the corresponding control point, u denotes the angle of the short edge 
with respect to the horizontal axis and x ^ specifies the chirality of the triangle. 

The pinwheel point set A forms a substitutive point set with substitution rule 
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where Rg = e*^ for G M, so multiplication by Rg corresponds to a rotation in the complex 
plane with the angle 9. A patch of the pinwheel tiling, with control points coloured according 
to the two chiralities of the triangular tiles, is shown in Figure [2j 
The infiation rule a gives rise to the substitution matrix 



M 



gitO _|_ giivr 2e**(2'p-7r) _|_ ^it{2if+^) 

1 + 2xy~^ + xy 

2y'^ + l/y^^ X + xy 



2 I 1 /„,-! ^ 1 ^„,-2 ; (2) 



where we set x = e^*"^* and y = e*'^*/^. This matrix was used in |12] to derive the circular 
symmetry of the autocorrelation. If we only wish to distinguish orientations modulo 90 
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Figure 2: A patch of the pinwheel tihng. The control point in A are coloured according to 
the two chiralities of the triangles. 

degrees, we can choose t = 4, so y = 1, and the matrix simplifies to 



This matrix can be used to derive information about the subsets of control points that belong 
to the same direction, in the sense that they belong to the same power of x. 

3 Hierarchies of subsets of different orientation 

Let us consider what happens when we start from an initial patch containing one triangle of 
each chirality with horizontal orientation, and perform n inflation steps. Denoting by fn{k) 
the number of triangles of each chirality at step n with orientation factor , where < < n, 
we find 




(3) 




(4) 



Multiplication by M gives the following recursion for the coefficients 



= 2/+(A;) + 3/-(A:-l), (5a) 
f-^^{k) = 2/+(A;-l) + 3/-(A:), (5b) 

for < /c < n + 1, where the appropriate initial conditions are /o^(0) = 1 and /^(A;) = for 
all fc < and k > n. 

Clearly, these numbers satisfy the relation fnik) = fn{^~ k) for all < k < n, as can 
be seen from the relation 

/O 1\ /2 3x\ /O 1\ f2x 3\ (2 2,x-^\ 

[l oj (,3 2x) [l oj = U V ^ ^3 2x-iJ- 

Hence the total number of triangles f^ik) = fn{k) + fn{k) obeys the symmetry fn{k) = 
fn{n — k) for all < A; < n. 

These observations have some interesting consequences. For instance, it can been shown 
that for even n 

/+(0) < /-(o) < /+(1) < /-(I) < . . . < /+(§ - 1) < /-(§ - 1) < /+(§) = /-(§ ) 
and for odd n 

/n+(0) < /„^(0) < /+(1) < /-(I) < . . . < /+(^ - 1) < /-(^ - 1) < /+(^) = /-(^) . 

It follows that fn{k) > fn{k — 1) for 1 < < [^] — 1, and the symmetry property 
fn{k) = /n(^ ~ k) then implies a corresponding inequality for large k. 



4 Diffraction of pinwheel point sets 

While it has been shown that the autocorrelation of the pinwheel point set is circularly 
symmetric [12], it is not known whether it is concentrated on sharp rings only, or whether 
there is also a continuous component. Direct numerical computations are of no help here, 
because any finite patch only contains a small number of independent directions, and is clearly 
not a good approximation of the limit structure. However, as shown in [Tl|2], the knowledge 
that the diff'raction is circularly symmetric can be exploited to derive an approximation of the 
diffraction based on a radial version of the Poisson summation formula; see \l\ for details. This 
results in an expansion in terms of Bessel functions. As shown in [Tl[2], this results in a radial 
distribution of the diffraction intensity that shows clear maxima for distances corresponding 
to those of a square lattice, so the radially discrete part resembles that of a powder diffraction 
of a square lattice. However, as the pinwheel structure contains a hierarchical sequence of 
scaled square lattices, the intensity distribution is not the same as for the square lattice, but 
deviates in a characteristic fashion. In addition, these numerical data provide an indication 
that there is also a continuous component. 

Here, we consider the slightly generalised case of diffraction from weighted pinwheel sets. 
Circular symmetry also holds for the autocorrelation of a weighted pinwheel control point 
set, with weights depending on the two chiralities of the triangular tiles. This follows from [9l 
Thm. 6.1]: the triangles can be regarded as two prototiles, a left-handed and a right-handed 
one, and each individual prototile then shows circular symmetry, too. 
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Figure 3: Part of the approximated radial intensity distribution for the pinwheel diffraction 
(o-l- = Q!_ = 1), compared with the exact diffraction for a square lattice powder (red bars). 



Denoting the sets of pinwheel control points for the two chiralities by A±, with Aj^ U /!_ 
A, we thus consider the weighted Dirac comb 



a+ 5a+ + a- 5a. ■= ^ a+ + ^ q_ 5^ , 



(7) 



where a± are arbitrary, in general complex, weights. In this paper, we restrict ourselves to 
values a± E {—1,0, 1}, with the case a+ = a- = 1 corresponding to the pinwheel diffraction 
considered in [H [2] . 

The approximations to the radial density distribution shown below are obtained by com- 
puting approximate autocorrelation coefficients r/(r) for distances r, obtained from finite ap- 
proximants of the pinwheel tiling, up to a suitable maximum distance. The approximation 
to the autocorrelation then has the form 
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where /i^ denotes the uniform measure on the unit circle, and 



V C 



5^ 



P,q,ie n} 



(8) 



(9) 



denotes the distance set for the pinwheel tiling. Here, the assumption of circular symmetry 
enters — the approximate autocorrelation coefficients obtained from a finite patch are used to 
define a circularly symmetric autocorrelation measure. From Equation ([8]), the approximate 
radial distribution of the diffraction intensity I{k) is obtained as the weighted sum of the 
Fourier transform of the uniform measure fir- As the latter is given by flr{k) = Jo(27r|/c|r), 
this results in finite sums of Bessel functions Jq which show strong oscillations; compare [IllS] 
for details. We note in passing that some values for frequencies of distances are known exactly, 
and they can in principle be computed from the inflation rule; see [13j . 



4.1 Systematics of peak intensities 

Figure [3] shows the numerical approximation for the diffraction of 5 a, which is the pinwheel 
point sets with weights a+ = a_ = 1. The corresponding result for a square lattice powder 
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Figure 4: Intensity ratios I{k)/I{l) of peaks for values of k such that k"^ is not divisible by 
5, as a function of k. On the left, different colours refer to different number of prime factors 
equal to 1 modulo 4 (either (green) , 1 (blue) or 2 (red) ) . Note that equal factors are counted 
here; for instance, the lowest k in the second group is k = \/l3, and in the third group it is 
k = \/l69 = 13. On the right, these data have been collapsed as described in the text. Lines 
are least square fits to functions of the form c/k with constant c. 
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Figure 5: Same as right part of Figure HI for values of k such that k'^ is divisible by 5 but not 
5^ (left) and such that k'^ is divisible by 5^ but not by 5^ (right). 

is also shown, with relative normalisation chosen to make the first peak match. Clearly, as 
observed in [Ti[2], the main peaks in the intensity are well reproduced, but there are systematic 
deviations in the peak intensities. The latter can qualitatively be explained as arising from 
distances in V of Equation ([9]) with denominators 5^ with i > 0. 

To get a better picture of the peak intensities, we consider separately the integrated 
intensities on rings at values of k according to the highest power of 5 that divides k'^ (which 
will be 0, 1 or 2 in the examples below). Figure H] shows the result for the case where k is 
not divisible by 5. On the left, you can clearly distinguish three groups of data, indicated by 
different colours, and for each group of data the decay of intensities with increasing k is well 
approximated by a function of the type c/k, where c is a constant. 

The three groups of data are distinguished by the number of prime factors of k'^ that are 
equal to 1 modulo 4. In fact, dividing the intensity ratios by 2^*, where s is the sum of powers 
of all prime factors of this type, collapses the data set onto a single curve, as shown in the 
right part of Figure [H This observation is in line with the similarity to the square lattice 
powder diffraction seen in [Tl[2]. 

The corresponding collapsed data for peaks at values of k which are multiples of 5 and 5^ 
are shown on the left and right of Figure El respectively. Again, the observed decay of the 
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Figure 6: Three approximations of the radial intensity distribution of the diffraction of the 
Dirac comb 5a+- 

intensities is inversely proportional to k. This indicates that our qualitative understanding 
of the peak intensities is correct, though it is not clear whether this also holds quantitatively 
for the diffraction of an infinite pinwheel tiling point set. 

4.2 Diffraction of weighted pinwheel point sets 

In order to get a better indication of the possible continuous contribution of the pinwheel 
diffraction, it is advantageous to consider weighted pinwheel point sets. As discussed above, 
the resulting autocorrelation is still circularly symmetric, and hence the same approach as 
outlined above can be applied to obtain approximations of the radial intensity distribution. 
Figure [6] shows the result for the case u = (a+ = 1 and a_ = 0), where only the half of 
the pinwheel point set is considered. The three spectra correspond to an increasing number of 
terms in the approximation, and increasing accuracy of the estimated correlation coefficients 
from larger patches of the pinwheel tiling constructed by successive inflation. The result very 
closely resembles that for the complete pinwheel point set, and may well be identical apart 
from the overall intensity normalisation due to the different density of scatterers. 

More interesting is the case of = 1 and q_ = — 1. Since this is a balanced case, the 
average scattering strength is zero, and there is no peak at k = 0. Therefore, a continuous 
background may be easier to spot in this case. Figure [7] shows three increasingly accurate 
approximations of the radial diffraction intensity for this case. 

There is not only an indication of the presence of a continuous background, similar to 
what was found in [2] , but also a number of new features appear to develop (in particular for 
small values of k), which may correspond to additional peaks for the infinite system. However, 
the approximation may not yet be good enough to draw firm conclusions. 
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Figure 7: Three approximations of the radial intensity distribution of the diffraction of the 
Dirac comb uj for the balanced case with a+ = 1 and a_ = — 1. 

5 Summary 

The diffraction of substitution tilings with continuous rotational symmetry still remains to 
be completely understood. In this paper, we presented some observations on combinatorial 
properties of the pinwheel tiling and on its diffraction measure, using the approximation 
introduced in [H [2]. From the numerical data, we corroborate the conclusions of [U |2] on 
the similarities between the pinwheel diffraction and a square lattice powder. Considering 
the balanced case of zero average scattering strength provides an indication that there might 
not only be a continuous contribution as conjectured in [2j, but that there might also be 
additional sharp rings with finite intensity. This observation warrants further investigation. 
It would also be interesting to compare this result with the diffraction of other tilings, such 
as the generalised pinwheel tilings introduced by Sadun [18J or the tipi tilings in 
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